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Abstract. In this paper, we are interested in sequences of q-tuple of N X N random matrices having 
a strong limiting distribution (i.e. given any non-commutative polynomial in the matrices and their 
conjugate transpose, its normalized trace and its norm converge). We start with such a sequence having 
this property, and we show that this property pertains if the q-tuple is enlarged with independent 
unitary Haar distributed random matrices. Besides, the limit of norms and traces in non-commutative 
polynomials in the enlarged family can be computed with reduced free product construction. This 
extends results of one author (C. M.) and of Haagerup and Thorbj0rnsen. We also show that a p-tuple 
of independent orthogonal and symplectic Haar matrices have a strong limiting distribution, extending 
a recent result of Schultz. We mention a couple of applications in random matrix and operator space 
theory. 

Dans cet article, nous nous interessons au q-tuple de matrices N X N qui ont une distribution limite 
forte (i.e. pour tout polyndme non-commutatif en les matrices et leurs adjoints, sa trace normalisee 
et sa norme convergent). Nous partons d'une telle suite de matrices aleatoires et montrons que cette 
propriete persiste si on rajoute au q-tuple des matrices indepandantes unitaires distribudes suivant la 
mesure de Haar. Par ailleurs, la limite des normes et des traces en des polyndmes non-commutatif s en 
la suite elargie peut etre calculi avec la construction du produit libre reduit. Ceci etend les resultats d 'un 
des auteurs (CM.) et de Haagerup et Thorbj0rnsen. Nous montrons aussi qu'un p-tuple de matrices 
indepandantes orthogonales et symplectiques ont une distribution limite forte, etendant par la-mime un 
resultat de Schultz. Nous passons aussi en revue quelques applications de notre resultat aux matrices 
aleatoires et a la theorie des espaces d'oporateur. 

1. Introduction and statement of the main results 

Following random matrix notation, we call GUE the Gaussian Unitary Ensemble, i.e. any sequence 
(X/v)tv>i of random variables where Xn is an N x N selfadjoint random matrix whose distribution is 
proportional to the measure exp ( — N/2Tr(A 2 )^dA, where dA denotes the Lebesgue measure on the set 
of N x N Hermitian matrices. We call a unitary Haar matrix of size N any random matrix distributed 
according to the Haar measure on the compact group of N by N unitary matrices. 

We recall for readers' convenience the following definitions from free probability theory (see |U|5T]). 
Definition 1.1. 

(1) A C*-probability space (A, .*, r, || ■ ||) consists of a united C* -algebra {A, .*, || • ||) endowed with 
a state t, i.e. a linear map r : A — > C satisfying t[1j\] = 1 and r[aa*] ^ for all a in A. In this 
paper, we always assume that t is a trace, i.e. that it satisfies r[ab] = r[ba] for every a,b in A. 
An element of A is called a (non commutative) random variable. A trace is said to be faithful 
if r[aa*] > whenever a =/= 0. If t is faithful, then for any a in A, 

(1.1) \\a\\= lim (r[(a*a) k ]). 

(2) Let A\, . . . ,Ak be *-subalgebras of A having the same unit as A. They are said to be free if for 
all ai £ Aj t (i = 1, . . . , k, ji G {1, . . . , k}) such that r[aj] = 0, one has 

r[a x • • • a k ] = 

as soon as j\ ^ ]2, ji ^ J3, ■ • • ,jk-l 7^ 3k- Collections of random variables are said to be free if 
the unital subalgebras they generate are free. 

(3) Let a = (a±, . . . ,a k ) be a k-tuple of random variables. The joint distribution of the family a is 
the linear form P >— > r[P(a, a*)] on the set of polynomials in 2k non commutative indetermi- 
nates. By convergence in distribution, for a sequence of families of variables (sln)n^i — 

l 
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(til \ . . . , ajj, )n^I in C* -algebras (An, •*, Tjv> || ■ ||), we mean the pointwise convergence of the 
map 

P^T N [P(a N ,a N )], 

and by strong convergence in distribution, we mean convergence in distribution, and point- 
wise convergence of the map 

P^\\P(a N ,a N )\\. 

(4) A family of non commutative random variables x = (x±, . . . , x p ) is called a free semicircular 
system when the non commutative random variables are free, self adjoint (xi — x* , i — 1. . . . ,p), 
and for all k in N and i = 1, . . . ,p, one has 

t[x*] = ( t k da(t), 



with do~(t) = — t 2 l|t|^2 dt the semicircle distribution. 

(5) A non commutative random variable u is called a Haar unitary when it is unitary (uu* — u*u — 
\a) and for all n in N, one has 

1 if n = 0, 
otherwise. 

In their seminal paper |15| . Haagerup and Thorbj0rnsen proved the following result. 

Theorem 1.2 ( [T5] The strong asymptotic freeness of independent GUE matrices). 
For any integer N ^ 1, let x[ N \ . . . ,Xp be N x N independent GUE matrices and let (xi, . . . ,x p ) 
be a free semicircular system in a C* -probability space with faithful state. Then, almost surely, for all 
polynomials P in p non commutative indeterminates, one has 

\\P(X[ N \...,XW)\\ — > llPOrx,...,^)!!, 

where \\ ■ \\ denotes the operator norm in the left hand side and the norm of the C* -algebra in the right 
hand side. 

This theorem is a very deep result in random matrix theory and had an important impact. Firstly it 
had significant applications to C*-algebra theory [151 122) . and more recently to quantum information 



theory (5J Ej • Secondly, it was generalized in many directions. Schultz (24] has shown that Theorem 1.2 
is true when the GUE matrices are replaced by matrices of the Gaussian Orthogonal Ensemble (GOE) 
or by matrices of the Gaussian Symplectic Ensemble (GSE). Capitaine and Donati-Martin [6J and, very 
recently, Anderson have shown the analogue for certain Wigner matrices. 

An other significant extension of Haagerup and Thorbj0rnsen's result was obtained by one author (C. 
M.) in [TS], where he showed that if in addition to independent GUE matrices, one also has an extra 
family of independent matrices with strong limiting distribution, the result still holds. 

Theorem 1.3 ( [TH] The strong asymptotic freeness of Xjy, Y^v). 
For any integer N ^ 1, we consider 

• a p-tuple Xjy of N x N independent GUE matrices, 

• a q-tuple Yjy of N x N matrices, possibly random but independent o/Xtv. 

The above random matrices live in the C* -probability space (Mjv(C), .*, tjv, | • ||), where tjv is the nor- 
malized trace on the set Mjv(C) of N x N matrices. In a C* -probability space (A, -*,T, || • ||) with faithful 
trace, we consider 

• a free semicircular system x of p variables, 

• a q-tuple y of non commutative random variables, free from x. 

If y is the strong limit in distribution o/Yjy, then (x, y) is the strong limit in distribution o/(Xat, Y^v). 

In other words, if we assume that almost surely, for all polynomials P in 2q non commutative indetermi- 
nates, one has 

(1-2) t n [P(Y n ,Y n )} — ► r[P(y,y*)] 5 

(1-3) UPCYw.YSrJll — > ||P(y,y*)||, 
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then, almost surely, for all polynomials P in p + 2q non commutative indeterminates, one has 
(1.4) ^[PCX^Y^Y^)] — > r[P(x,y,y*)], 

Jy — too 

(1-5) ||P(X N ,Y 7V ,Y^)|| — ► ||P(x,y,y*)||. 



The convergence in distribution, stated in ( 1.4 1, is the content of Voiculescu's asymptotic freeness theorem. 
We refer to [U Theorem 5.4.10] for the original statement and for a proof. An alternative way to state 
(1.5 1 is the following interversion of limits: for any matrix Hn = _P(Xjv , Yjv, Y|^), where P is a fixed 
polynomial, if we denote h — P(x, y, y*), then by the definition of the norm in terms of the state (1.1), 

lim lim (T N [(H* N H N ) k ]Y k = lim (r[(h*h) k ] V* . 

JV->oo k-^oo \ J fe->oo V L J 

It is natural to wonder whether, instead of GUE matrices, the same property holds for unitary Haar 
matrices. The main result of this paper is the following theorem. 

Theorem 1.4 (The strong asymptotic freeness of u\ N \ . . . , Up , Yjv). 

For any integer N J? 1, we consider 

• a p-tuple Ujv of N x N independent unitary Haar matrices, 

• a q-tuple Yjy of N x N matrices, possibly random but independent o/Ujv- 

In a C* -probability space (A, .*, r, |j • ||) with faithful trace, we consider 

• a p-tuple u of free Haar unitaries, 

• a q-tuple y o/ non commutative random variables, free from u. 

7/y «s i/ie strong limit in distribution o/Yjv, £/ien (u, y) is the strong limit in distribution o/(Utv, ~Yn)- 

In order to solve this problem, it looks at first sight natural to attempt to mimic the proof of [T5] and 
write a Master equation in the case of unitary matrices. However, even though such an identity can be 
obtained for unitary matrices, it is very difficult to manipulate it in the spirit of [15 in order to obtain 
the desired norm convergence. Part of the problem is that the unitary matrices are not selfadjoint, unlike 
the GUE matrices considered in [15], and in this context the linearization trick and the identities do not 
seem to fit well together. In order to bypass this problem, in this paper, we take a completely different 
route by building on Theorem |1.3| and using a series of folklore facts of classical probability and random 
matrix theory. 

Our method applies to prove the strong convergence in distribution of Haar matrices on the orthogonal 



and the symplectic groups by building on the result of Schultz [24] , which is the analogue of Theorem 1.2 
for GOE or GSE matrices instead of GUE matrices. The analogue of Theorem |1.3| does not exist yet. If 
one shows that the estimates of matrix valued Stieltjes transforms in |19| can always be performed with 
the additional terms in the estimate of [24], then, following the lines of this paper, one gets Theorem 



1.3 for Haar matrices on the orthogonal and the symplectic groups, instead of the unitary group only. 
Therefore, in the following Theorem, we stick to proving the strong convergence of independent unitary, 
orthogonal or symplectic Haar matrices, without "constant" matrices Y: 

Theorem 1.5 (The strong asymptotic freeness of independent Haar matrices). 

For any integer N ^ 1, let u[ , . . . , Up be a family of independent Haar matrices of one of the three 
classical groups. Let u\, . . . , u p be free Haar unitaries in a C* -probability space with faithful state. Then, 
almost surely, for all polynomials P in 2p non commutative indeterminates, one has 



\p(ui N \...,uj ) N \u^*,...,uj 1 N >)\\-^ \\P( U1 ,. 



N- 

where \\ ■ \\ denotes the operator norm in the left hand side and the C* -algebra in the right hand side. 

Our paper is organized as follows. Section [2] consists of applications of the results stated above. Among 
other examples, we show that the limit of complicated random matrix models involving unitary random 
matrices have norms that converge towards (or are bounded by) values predicted by the theory of free 
probability. Sections [3] and [4] provide the proofs of Theorem 1 1 . 4| and Theorem |1.5| respectively. Section [5] 
is dedicated to the proof of Corollary |2.2[ stated in the next section. 
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2. Applications 
Our main result has the potential for many applications. 
2.1. The spectrum of Hermitian random matrices. 

2.1.1. Generalities on the strong convergence in distribution. We first recall for convenience some facts 
about the strong convergence in distribution, mainly an equivalent formulation. Given a self-adjoint 
variable h in a C* -probability space (A, *,t, || • ||), its spectral distribution //^ is the unique probability 
measure that satisfies r[h k ] = J t k d^i(t) for any k ^ 1. This measure has compact support included in 
[ — \\h\\, \\h\\\ . For any continuous map / : R to C, the variable f(h) is given by functional calculus, and 
coincides with the limit of (Pn(h)) n>1 in A, where (P n ) n ^i is an Y Weierstrass's approximation of / by 

polynomials. Given a (non self-adjoint) variable x in A, we set the self-adjoint variables Jfx = (^ x+ 2 ^ 
and 3x = {^^-^j , so that x = + iQx. 

Proposition 2.1 (The strong convergence in distribution of self adjoint random variables). 

Let xat = (x[ N \ . . . , Xp N ^) and x = (xi,...,x p ) be p-tuples of variables in C* -probability spaces, 

(An, ■* ,t~n, II • II) and (A, .*, r, || • ||), with faithful states. Then, the following assertions are equivalent. 

(1) xat converges strongly in distribution to x, 

(2) for any continuous map fi,gi : M — > C, i = l,...,p, the family of variables 
(fi(^x[ N ^),gi(^x[ N ^), . . . , fp(^x p N ^), gp^x^)) converges strongly in distribution to 
(/i(SRa;i),ffi(9tei), . . . , fp{9tx p ),g p (5Xx p )), 

(3) for any self-adjoint variable hpj — P(xjv), where P is a fixed polynomial, fih N converges in weak-* 
topology to \i h where h = P(x). Weak-* topology means relatively to continuous functions on C. 
Moreover, the support of Hh N converges in Hausdorff distance to the support of ji^, that is: for 
any e > 0, there exists No such that for any N Nq, 

(2.1) Supp( n hN ) C Supp( [L h ) +(-e,e). 

The symbol Supp means the support of the measure. 
In particular, the strong convergence in distribution of a single self-adjoint variable is its convergence in 
distribution together with the Hausdorff convergence of its spectrum. 

Proof. Assuming (1), the assertion (2) is obtained by Weierstrass's approximation of the functions /$ and 
<7i by polynomials in p complex variables on the centered ball of radius sup^g ||a;jv||. The converse is 
obvious. 



Assuming (1), let us show (3). By Weierstrass's approximation, converges strongly in distribution to 
h. The convergence in distribution of h?q to h implies the weak-* convergence of ^ N to /i^. For any e > 0, 
let f e be a continuous map which takes the value 1 on the complementary of Supp ((J-h) + (~ £ : £ ) an d 
on Supp (fj, h ). Then, ||/ £ (x w )|| — > ll/ £ ( x )ll = ^k(J / £ (x) fe d/i /l )s = 0. Hence, the support of fj, hN 

N— >oo 

is a subset of Supp + £, e) for N large enough, as otherwise one could have ||/ e (xjv)|| = 1 eventually. 

Assuming (3), let us show (1). Let P be a polynomial in p variables and their conjugate. Denote — 
P(x N ,x* N ) and m = P(x,x*). Then, T N [m N ] - r[m] = t n [Q(x n , x* n )] - r[Q(x, x*)] +i(r N [R(x N , x* N )] - 
t[R(x, x*)]) where Q = \(P + P*) and R = ^(R — R*) gives Hermitian variables. By the assertion (3) 
and since the matrices are uniformly bounded in operator norm, we get the convergence in moments of 
the spectral distribution of Q(xn,x* n ) and R(xn,x* n ). Hence, we get the convergence in distribution of 
xn to x. Then, the convergence holds in weak-* topology since /i/j has bounded support. Furthermore, 

||mAr|| 2 = ||m^mjv|| = maxSupp ( £t m * mjv ) — > maxSupp ( [i m » m ) = \\m*m\\ = \\m\\ 2 . 

™ N— >ca 

□ 

2.1.2. The spectra of the sum and product of unitary invariant matrices. The following is a consequence 
of our main result: 

Corollary 2.2. Let An,Bn be two N x N independent Hermitian random matrices. Assume that: 
(1) the law of one of the matrices is invariant under unitary conjugacy, 
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(2) almost surely, the empirical eigenvalue distribution of An (respectively Bn) converges to a com- 
pactly supported probability measure fj, (respectively v), 

(3) almost surely, for any neighborhood of the support of [i (respectively v), for N large enough, the 
eigenvalues of An (respectively Bn) belong to the respective neighborhood. 

Then, one has 

• almost surely, for N large enough, the eigenvalues of An + Bn belong to a small neighborhood of 
the support of \i ffl v, where EH denotes the free additive convolution (see |211 Lecture 12]J. 

• if moreover Bn is nonnegative, then the eigenvalues of (B/v) 1 / 2 An(Bn) 1 ^ 2 belong to a small 
neighborhood of the support of fiMv, where Kl denotes the free multiplicative convolution (see |21l 
Lecture 14]/ 



Corollary 2^ is proved in Section [5] It can be applied in the following situation. Let An be an N x N 



Herniitian random matrix whose law is invariant under unitary conjugacy. Assume that, almost surely 
the empirical eigenvalue distribution of An converges to a compactly supported probability measure [i 
and its eigenvalues belong to the support of [i for TV large enough. Let n^r be the matrix of the projection 
on first pn coordinates, Hn = diag (1 PJV , 0n- Pn ), where pn ~ tN, t G (0, 1). We consider the empirical 
eigenvalue distribution fiN of the Herniitian random matrix 

Then, it follows from a Theorem of Voiculescu [M] (see also (JJ) that almost surely [In converges weakly 
to the probability measure ^ft = [i S3 [(1 — t)So + t8i]. This distribution is important in free probability 
theory because of its close relationship to the free additive convolution semigroup (see |211 Exercise 
14.21]). Besides, the empirical eigenvalue distribution /ijv was proved to be a dctcrminantal point process 
obtained as the push forward of a uniform measure in a Gelfand-Cetlin cone [10]. Very recently, it was 
proved by Metcalfe [20 that the eigenvalues satisfy universality property inside the bulk of the spectrum. 
Our result complement his, by showing that almost surely, for N large enough there is no eigenvalue 
outside of any neighborhood of the spectrum of [ft' . 

2.2. Questions from operator space theory. We present some examples of norms of large matrices 
we can compute by Theorem |1.4| as the norm of the limiting variables have been computed by other 
authors. 



2.2.1. The norm of the sum of unitary Haar matrices. The following question was raised by Gilles Pisier 
to one author (B.C.) ten years ago: let U[ N 
matrices, p ^ 2. Is it true that almost surely: 



to one author (B.C.) ten years ago: let u[ N \ . . . , Up N ^ be N x N independent unitary Haar random 



v 



(2-2) ,77* 2V^T. 



JV- 



This question is very natural from the operator space theory point of view |22l Chapter 20] , and was still 
open before this paper. Haagerup and Thorbj0rnsen's theorem |15| have proved that the the convergence 
( pT2| is true when u[ N \ U^ N) are certain sequence of independent large unitary matrices (non Haar 



distributed). Our main theorem implies that (2.2) is true almost surely when they are i.i.d. unitary Haar 
matrices. Indeed, 2^/p — 1 is the norm of the sum of p free Haar unitaries by a result of Akemann and 
Ostrand [Ij: they have proved that if Ui are the generators of the free group von Neumann algebra, then 

p , p . 

(2.3) J) E a *^ = mm |2t + £j (V* 2 + M 2 - t) j. 



And if a\ — . . . = a p = 1 they prove that the minimizer of the right hand side is 2^/p — 1. 



By Theorem 1.5 and (2.3 1, we get that, for independent Haar matrices U± , . . . , Up N ^ on the orthogonal 



unitary or symplectic group, almost surely one has 

II £ aiU * N) II ^ i 2t + £ (v^KF - *) }, 



which is a generalization of (2.2) 
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2.2.2. The sum of Haar matrices along with their conjugate. In the same vein, by a result of Kesten |17| . 
the norm of the sum of p free Haar unitaries and of their conjugate equals 2^/2p — 1. Hence, we get from 
our result that almost surely one has 



(AO 



u. 



(AT)* 



2yj2 P -\. 



Remark that this result is not true for random unitary matrices distributed according to the uniform 
measure on the set of permutation matrix. Indeed, in that case 2p is always an eigenvalue of the matrix 
since N ^ + *) is the adjacency matrix of a 2p-regular graph. The convergence of the second 

largest eigenvalue to 2^/2p — 1, known as Alon's conjecture [2], has been proved recently by Friedman 

ED. 



2.2.3. The sum of Haar matrices, matrix valued case. Lehner |18| has proved that for m, . . . ,u p free Haar 
unitaries and a ,ai, . . . ,a p Hermitian k by k matrices 



(2.4) 









Ui 


= inf 






6>0 





where the infimum is over all positive definite invertible k by k matrices b. Recall that from Theorem 1.5 
we can deduce the following corollary (see |191 Proposition 7.3] for a proof). 

Corollary 2.3. Let U be a family of independent Haar matrices one of the three classical groups. Let u 
be a family of free Haar unitaries. Let k 1 be an integer. Then, for any polynomial P with coefficients 
in Mfc(C), almost surely one has 



\\P{V N) \J$ 



N J 



N- 



\P(n,u*)\\, 



and in the right hand side for 



where \\ ■ \\ stands in the left hand side for the operator norm in M^nI 
the C* -algebra norm in Mfc(_4). 

We then deduce that the norm of block matrices of the form oq <8) 1 + Ya=i a i ® where ao, 

are Hermitian matrices, converges almost surely to the quantity (2.4 1 computed by Lehner. 



. ax 



2.2 A. Application of Fell's absorption principle. For an other application of Corollary |2.3[ recall Fell's 
absorption principle |22l Proposition 8.1]: for any k by k unitary matrices a\, . . . ,a p and m, . . . ,u p a free 
Haar unitaries, one has 

, p 



^ eg ® Ui = y^ttJ = 2^Jp - 1. 



By Corollary 2.3 we get for any k x k unitary matrices a\, . . . ,a p , almost surely one has 



53 a* <g> U, 



(AO 



N- 



which solves a question of Pisier in [221 Chapter 20] . 
2.3. Estimates on the norm of random matrices. 

2.3.1. Haagerup's inequalities. Let u = (tti, . . . ,u p ) be free Haar unitaries in a C* -probability space 
(A, .*, r, || • ||) with faithful state. For any integer d ^ 1, we denote by Wd the set of reduced *-monomials 
in p indeterminates x = [x\, . . . , x p ) of length d: 



W d 



P : 



3l 



3d 



3i 3d, e, G {1,*} Vj = l,...,d }. 



In 1979, Haagerup [14] has shown that one has 



(2.5) 



J2a n Pn{u) <(d + l)||a|| 2 , 



for any sequence (P n ) n ^i of elements in Wd and sequence a 
€ 2 -norm is denoted by 



(cKn)n^i of complex numbers whose 



"2 = 



it, 
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This result, known as Haagerup's inequality, has many applications (for example, estimates of return 
probabilities of random walks on free groups) and has been generalized in many ways. For instance, 



Buchholz has generalized (2.51 in an estimate of 5Z n >i a « ® x ni where the a n are now k x k matrices. 



Let Uat be a family of p independent N x N unitary Haar matrices. As a byproduct of our main result, 
we get that 

lim ||y>„P„(U w )|| ^ (d+l)||a|| 2 , 



N- 



where for any n 1, the polynomial P n is in Wd- 

2.3.2. Kemp and Speicher's inequality. Kemp and Speicher |16| have generalized Haagerup's inequality 
for 7?.-diagonal elements in the so-called holomorphic case (with polynomials in the variables, but not 



their adjoint). Theorem 1.4 established, the consequence for random matrices sounds relevant since it 



allows to consider combinations of Haar and deterministic matrices, and then get a bound for its operator 
norm. The result of [16] we state below has been generalized by de la Salle (9] in the case where the non 
commutative random variables have matrix coefficients. This situation could be interesting for practical 
applications, where block random matrices are sometimes considered (see |25| for applications of random 
matrices in telecommunication). Nevertheless, we only consider the scalar version for simplicity. 

Recall that a non commutative random variable a is called an 7?.-diagonal element if it can be writ- 
ten a = uy, for u a Haar unitary free from y (see |21|). Let a = (ai, . . . , a p ) be a family of free, identically 
distributed 7?.-diagonal elements in a C* -probability space (A, *,t, || ■ ||). We denote by the set of 
reduced monomials of length d in variables x (and not its conjugate), i.e. 



°0d 



hi 1 -- - ¥=jd }■ 



Kemp and Speicher have shown the following, where the interesting fact is that the constant (d + 1) is 
replace by a constant of order y/d+l: for any sequence (P n ) n ^i of elements of and any sequence 
i = (a n )n^i, one has 



(2.6) 



^2 a nPn(&) «S eVd+ 1 ^ a nPn{a.) 



where || • H2 denotes the L 2 -norm in A, given by \\x\\2 = t^*^] 1 / 2 for any a in A. In part icula r, if a = u 
is a family of free unitaries (i.e. y = 1) then we get || J2 n >i a nPn(u) 



an improvement of (2.5 1 without the generalization on 7?.-diagonal elements 



a\\2, so that (2.6 1 is already 



N 



Vp ) be families of N x N independent unitary 
Yp"') be a family of N x N deterministic Hermitian matrices. As- 

(N) 



(N) 



Now let V N = (Ul 
Haar matrices and Yjy 

sume that for any j — 1, . . . ,p, the empirical spectral distribution of Yj " ' converges weakly to a measure 

fi (that does not depend on j) and that for N large enough, the eigenvalues of Y^ belong to a small 
neighborhood of the support of [i. We set for any j — 1, . . . ,p the random matrix 



A {N) = U (N) Y {N) V {N) * 



From Theorem 



1.4 



and [19, Corollary 2.1], we can deduce that almost surely the family A^r = (A± 



(AT) 



..4 



converge s str ongly in law to a family a of free 7?.-diagonal elements, identically distributed. Hence, in- 
equality (2.6 1 gives: for any polynomials P n in , n ^ 1, 



lim || V a n P n (A[ N \ . . . ,4^)|| < eVd+l\\ V a n P n (a) 



>1 



3. Proof of Theorem 11.41 

We consider a unitary Haar matrix Un, independent of a family of matrices Yjv, having almost surely 
a strong limit in distribution. We show that almost surely {Un,~Yn) has almost surely a strong limit in 
distribution. As it is known that (t/jy, Yjv) converges in distribution [4, Theorem 5.4.10], the only thing 
we have to show is the convergence of norms. This will show Theorem |1.4| by recurrence on the number 
of Haar matrices. Moreover, the problem can be simplified in the following way (see [191 Section 3]): 
• one can reason conditionally, and then assume that the matrices of are deterministic, 



s 
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• one may assume that the matrices of Yjv are Hermitian by considering their Hermitian and 
anti-Hermitian parts, 

• it is sufficient to prove that for any polynomial P, almost surely the norm of ||P(Ujv, Ujv, Yjv) I 
converges, rather than "almost surely, for any polynomial". 

The keystone of the proof is the use of a classical coupling of real random variables, namely the inverse 



transform sampling method, for Hermitian matrices (Lemma 3.1 below). It allows us to get the strong 
convergence of (C/jv, Yjv) from the strong convergence of (Xn, Yjv), where Xn is a GUE matrix indepen- 
dent of Yjv, for which we know the strong convergence by Theorem |1. 3 1 For that purpose, we will first go 
through an intermediate problem. We use the coupling to prove in Lemma |3. 2 1 the strong convergence of 
(Mjv , Yjv ) , where Mjv is the unitary invariant random matrix whose spectrum is { , . . . , N ^ , ^ } ■ From 
this, we deduce that the strong convergence holds for (Zn, Yjv), where Zn is any unitary invariant ran- 
dom matrix, independent of Yjv, whose spectrum is {~/N(jj), ■ ■ ■ , Jn(^n~^), 1n(jj) } for 7jv : [0, 1] — >• C 
a random map converging uniformly to a continuous map. We finally use the coupling method again, to 
remark that a unitary Haar matrix could be written as above. 

3.1. Preliminaries. Let a be a self-adjoint element in a C*-algebra, that is a* — a. Denote by \i a 
its spectral distribution, i.e. fi a is the unique probability measure on K such that for any k ^ 1, 
/ < fc d/x (i) = r[a k ]. This measure has compact support included in [ — ||a|| , ||a||] . Denote by F a its 
cumulative function, satisfying, for all t in M: 

(3.1) F„(t):=M(-]oo,t]). 
We set the generalized inverse of F a : for any s in ]0, 1] 

(3.2) F-\s) = inf {t G [-tt.tt] | F a {t) > s}. 

By the inverse method for random variables |12| Chapter two] , we get the following lemma. 

Lemma 3.1 (The coupling of self-adjoint variables and Hermitian matrices by cumulative functions). 

(1) Let a, b be two self-adjoint non commutative variables. Denote the self-adjoint variables b = 
F^ o F a (a) given by functional calculus. If /i a have no discrete part (i.e. /i a ({t}) = for any t 
in R), then b has the same distribution as b. 

(2) Let An and B>n be two Hermitian matrices (living in the C* -probability space 
(Mjv(C),.*,tjv,|| • ID/ Write the matrices A N = V An A An V^ n and B N = V Bn A Bn V^ n , 
with Va n , Vb n unitary matrices, such that the entries of Aa n and Ab n are non decreasing along 
the diagonal. Assume that diagonal entries of Aa n are distinct. We set the matrix 

, s /I N —1 N\ . 

(3-3) M N := V An diag — — , ^) V* An . 



~N ' N N 

Then Mjv = F An {A n ) and F^(M N ) = V An A Bn VX n - 

3.2. Step 1: from the GUE to an intermediate model. 

Lemma 3.2 (The strong asymptotic freeness of Mm, Yjv). 
Define the random matrix 

/ 1 N - 1 N\ 
(3-4) Mjv = Fjv diag (-,... -J V N , 

where Vn is a unitary Haar matrix, independent of Yjv- Then, almost surely (Mjv,Yjv) converges 
strongly in distribution to (m,y), where y is the strong limit of Yjy, free from a self adjoint variable m 
whose spectral distribution is the uniform measure on [0, 1] . 

Proof. Let Xn be a GUE matrix independent from Yjv, such that Xn = VjvAjvVjv, where Ajv is a 
diagonal matrix, independent of Vn, with non decreasing entries along the diagonal (we recall a proof 
of that decomposition in Proposition 6.1 Section [6j. Let i be a semicircular variable free from the 



strong limit y of Yjv- Let Fx N and F x be the cumulative functions of the spectral measures of Xn and 



x respectively. By the coupling of Lemma 3.1 we get that m = F x (x) has the expected distribution 
and, since the eigenvalues of a GUE matrix are almost surely distinct, we get that almost surely Mjv = 
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Fx n (Xn)- Then, for any polynomial P, almost surely 



(3.5) 



\P(m,y)\\-\\P(M N ,Y N )\\ 



\P(F x (x),y)\\-\\P(F x (X N ),Y N )\\ 
- \\P(F X (X N ),Y N ) - P(F Xn (X n ),Y n )\\. 



The first term in the right hand side of (3.5) tends to zero almost surely by the strong convergence in 
distribution of (Xx,Yx) to (x, y) (Theorem 1.3 1 and Proposition 2.1 since F x is continuous. For the 



second term, recall first that the convergence in distribution of Xx to x implies the pointwise convergence 
Fx N to F x (at any point of continuity of F x , and so on E). By Dini's theorem (231 Problem 127 chapter 
II], Fx N convergence actually uniformly to F x . Hence, since the matrices Xx, Y x are uniformly bounded 
in operator norm and by local uniform continuity of P, the second term converges also to zero. □ 

3.3. Step 2: from the reference model to other unitary invariant models. 

Lemma 3.3 (The strong asymptotic freeness of Zx, Yx)- 
Consider an N by N random matrix Zx of the form 



(3.6) 



Z 



N 



7A r(M 



N 



) = V N diag ^ 7JV (1),..., 7W (^_I), 7JV ^^ Vj 



where Mx is the random matrix of Lemma 3.2 and "fx ■ [0,1] — > C is a random map, independent of 
Yx- Assume that almost surely "fx converges uniformly to a continuous map 7 : [0, 1] — > C, that is 

||t — TjvIIoo := sup |t(*) -7iv(*)| 



te[o,i] 



JV- 



0. 



We set the self adjoint variable z = 7(771) given by functional calculus, with m as in Lemma 3.2 (it is well 
defined since \\m\\ ^ 1). Then, almost surely, (Zx,Yx) converges strongly to (z.y). 

Proof. For any polynomial P, one has 



\P(z,z*,y)\\-\\P(Z N ,Z N Y N ) 



(3.7) 

(3.8) 
(3.9) 

where 7 denotes the complex conjugacy of 7. The first term of the right hand side of (13 .Th tends to 



P( 7 (m),7(m),y)|| - \\P(y(M N )^(M N ),Y N ) 
+ \\P( 1 (M N )n(M N ) 1 Y N ) - P(in(M n ),1n(M n ),Yx)\ 



zero by Lemma |3.2| Proposition |2.1| and the continuity of 7. By the uniform convergence of 7^, the 
continuity of polynomials, and the fact that the matrices Mx,Yx are uniformly bounded in operator 
norm, the second term vanishes at infinity. □ 

3.4. Step 3: application to Haar matrices. Now, let Ux be a unitary Haar matrix, independent 
of Yx- By the spectral theorem (see Proposition 6.1 in Section [6]), we can write Ux = Vn&nV n , 
where the entries of Ax = diag [e 2 ^ 101 



diagonal. Define as above Mx = Vx diag (-^, . . . , N N ' 1 , ^) V N , where Vn is the unitary matrix in the 



decomposition of Ux- Denote by Fjj n the cumulative function of diag \y l 
coupling of Lemma |3.1| that 

(3.10) 



!jv ) have non decreasing argument in [0,27r[ along the 

ary matrix in the 

). We get by the 



(N) 



9^ 



U N =exp(27riF c 7j(M J v)). 



By Lemma 3.3 to get the strong convergence of (Ux,Yx) it remains to prove that almost surely 
jx '■ t — > exp (2mFu^ (i)) converges uniformly to 7 : t — ¥ exp (2ivit). 

From the convergence of Ux to a Haar unitary u, we get that almost surely Fjj n converges to F u . 
Let t in [0, 1[. Almost surely, for any < a < 1 — t, there exists No ^ 1 such that for any N ^ No, 
Fu N {t + ct) ^ t + §. The points (Fy^(t),t) and (t + a,Fu N {t + a)) belong to the graph of Fu N , with 
vertical segments on points of discontinuity. Hence, since Fu N is non decreasing we get F^(t) ^ t + a. 
Hence, for any ^ t < 1, we get limsup^^^ F^ 1 (t) ^ t. With a symmetric reasoning we get 



liminfjv^ooF^i) ^ t. Now, remark that F^(0) = B x ^ and F^(l) 



'N 



^ 1. Hence, F, 



U N 



converges pointwise to the identity map on [0, 1]. By Dini's theorem, it converges uniformly. Hence 7jv 
converges uniformly to 7. 
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4. Proof of Theorem 11.51 

The proof of |1.5| is obtained by changing the words unitary, Hermitian and GUE into orthonormal, sym- 
metric and GOE. respectively symplcctic, self dual and GSE, and by taking Yn a family of independent 



orthogonal, respectively symplectic, matrices. Instead of Theorem 1.3 we use the main result of [24]. In 
the symplectic case, we also have to consider matrices of even size. 



5. Proof of Corollary 12.2 



First recall the following consequence of [19, Corollary 2.1]. 

Lemma 5.1. Let d[ N ^ and D2 be two diagonal matrices having a strong limit in distribution separately. 



Then, there exists diagonal matrices D\ 



(N) 



J 1 



with the same eigenvalues as D[ N ^ and D 



respectively, such that (D\ , D\ ) converges strongly in distribution. 



Let An and Bn be as in Corollary |2.2| Without loss of generality, we can assume that the laws of and 
Bn are invariant under unitary conjugacy. Let Aa n and Ab n be diagonal matrices of eigenvalues of An 
and Bn respectively. By Proposition |2.1| the assumptions on An and Bn mean their strong convergence 
in distrib ution separately, and so the strong convergence of A,4 N and A b n separately. With the notations 

Let (Un, 



of Lemma 
of (A An , 



5.1 



consider and Ag, 



A Bn ). Then (A N ,B N ) and (U N A, 



, Vjv) be independent unitary Haar matrices, independent 
N V N ) are pairs of random matrices with the 
we get the almost sure strong convergence 



r U N ,V N A B 



1.4 



same probability law (see Proposition 6.1). By Theorem 
of (Un 7 Vn 7 Aa n , Ab n ), and then of (UnAa n U n ,VnAb n V n )- Hence, we obtain that (An,Bn) has a 
strong limit in distribution (a, b). The spectral distribution of a is fx, the one of b is v, and a and b are 
free. The strong convergence implies the convergence of the spectrum of An + Bn to the support of jj, EB v 



(which is the spectral distribution of a + b) by Proposition 2.1 We then get the first point of Corollary 2.2 



2.2 



We get the second point of Corollary 
stated after Corollary |2.2| follows by taking 11^ 
and remarking that 11^ = II jv- 



1 li 

with the same reasoning on (A^^A^). The application 
Bn, which satisfies the assumptions since t € (0, 1), 



6. Appendix: The spectral theorem for unitary invariant random matrices 

This result seems to be folklore in the literature of Random Matrix Theory, but we are not able to find 
an exact reference, so we include a proof for the convenience of the readers. 

Proposition 6.1 (The spectral theorem for unitary invariant random matrices). Let Mn be an N x N 

Hermitian or unitary random matrix whose distribution is invariant under conjugacy by unitary matrices. 
Then, Mn can be written Mn — VnAnV n almost surely, where 

• Vn is distributed according to the Haar measure on the unitary group, 

• An is the diagonal matrix of the eigenvalues of Mn, arranged in increasing order if Mn Her- 
mitian, and in increasing order with respect to the set of arguments in [— 7r,7r[ if Mn is unitary, 

• Vn and An are independent. 

We actually use the proposition only for unitary Haar and GUE matrices, which are two cases where 
almost surely the eigenvalues are distinct. The fact that multiplicities of eigenvalues is almost surely one 
brings slight conceptual simplifications in the proof, but nevertheless do not change the result. Hence, 
we choose to state the proposition without any restriction on the multiplicity of the matrices. 

Proof. By reasoning conditionally, one can always assume that the multiplicities of the eigenvalues of Mn 
are almost surely constant. We denote by (Ni, . . . , Nk) the sequence of multiplicities when the eigenval- 
ues are considered in the natural order in R or in increasing order with respect to their argument in [—it, 7r[. 

Since Mn is normal, it can be written Mn — VnAnVn, where Vn is a random unitary matrix and 
An is as announced. The choice of Vjv can be made in a measurable way, see for instance |11[ Section 
5.3], with minor modifications. 

Let (iti, . . . ,uk) be a family of independent random matrices, independent of (An,Vn) and such that 
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for any k = 1, . . . , K, the matrix is distributed according to the Haar measure on U(N/.), the group 
of Nk x Nk unitary matrices. We set 

V N = V N diag (ui, . . . , u K ), 

and claim that the law of Vn depends only on the law of Mn, not in the choice of the random matrix 
Vn- Indeed, let Mn — VnAnVn be an other decomposition, where Vn is a unitary random matrix, inde- 
pendent of (iti, . . . , uk )• The multiplicities of the eigenvalues being N±, . . . , Nk, there exists (v%, . . . , vk ) 
in U(Ni) x • • ■ x U(Nk), independent of . . . , uk), such that Vn = Vn diag (yx, ■ . . , vk)- Hence, we 
get Vn diag (ui, . . . , Uk) — Vn diag (villi, . . . , vkUk), which is equal in law to Vn- This proves the claim. 

Let W N be an N x N unitary matrix. Then W N M N W N = (WnVn)An(W n Vn)*- By the above, 
since Mn and WnMnW n are equal in law, then Vn and WjvViv are also equal in law. Hence Vn is Haar 
distributed in U{N). 

It remains to show the independence between Vn and An- Let / : U(N) — > C and g : Mjv(C) — > C two 
bounded measurable functions such that g depends only on the eigenvalues of its entries. Then one as 
E[/(V;v)ff(Ajv)] = E[f(V N )g{M N )]- Let W N be Haar distributed in U(N), independent of (V N ,A N ). 
Then by the invariance under unitary conjugacy of the law of Mn, one has 

E[f(V N )g(A N )] = E[f(W N VN)g(WNM N W N )] 

= E[f(WNVN)g(A N )] 

= E[E[f(W N VN)\VN,A N ]g{AN)] 

= E[f(W N )]E[g(A N )] =E[f(VN)]E[g(A N )]. 

□ 
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